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Abstract 

We present the distribution of virial masses for nearby galaxy clusters, as 
obtained from a data-set of 75 clusters, each having at least 20 galaxy mem- 
bers with measured redshifts within 1.5 Mpc. After having accounted 
for problems of incompleteness of the data-set, we fitted a power- law to the 
cluster mass distribution. 
Subject headings: galaxies: clustering 



1 Introduction 

The distribution functions of observational cluster quantities, such as velocity 
dispersions, radii, masses, luminosities and X-ray temperatures, can provide 
strong constraints both on cosmological scenarios and on the internal dy- 
namics of these systems. Theoretical as well as observational estimates of 
these distribution functions are presently being debated. 

The theoretical mass distribution expected for groups and clusters in the 
hierarchical clustering scenario has been derived by several authors (see, e.g.. 
Press & Schechter 1974; Cavaliere, Colafrancesco, & Scaramella 1991; Blan- 
chard, Valls-Gabaud, & Mamon 1992). The distribution functions for mass. 
X-ray temperature and velocity dispersion of rich clusters have been shown 
to provide important diagnostics for cosmological models (see, e.g.. Edge et 
al. 1990, Henry & Arnaud 1991, Lilje 1992). In particular, the distribution 
function of cluster X-ray temperatures has been used to constrain the mass 
fluctuation spectrum (see, e.g., Henry & Arnaud 1991). Bahcall (1979) es- 
timated the cluster luminosity function, and first offered the possibility of 
evaluating the cluster mass function by adopting a constant mass-to-light 
ratio. At present both optical and X-ray data contribute to the description 
of the cluster mass function (Bahcall & Cen 1992). Recently, Pisani et al. 



(1992) obtained the mass distribution function of groups of galaxies directly 
from their dynamics, while Ashman, Persic, & Salucci (1993) described a 
distribution function of galaxy masses. 

Quite a large number of redshifts for cluster galaxies is now available, 
so that many cluster masses can be measured directly from dynamics. We 
collected from the hterature 75 clusters, each having at least 20 members 
with measured redshift within 1.5 Mpc (we use Hq — 100 h km Mpc~^ 
throughout). In Table 1 we list the clusters considered, the number N of 
cluster members with available redshift, and the richness classes R mainly 
taken from Abell, Corwin, & Olowin (1989; AGO hereafter). 

In order to reasonably reduce evolutionary effects, we considered only 
clusters with mean redshift z < 0.15. In Biviano et al. (1992), Girardi et al. 
(1993a), and Girardi et al. (1993b) we describe the data-base and we give 
the relevant references. In Girardi et al. (1993a) we detail the criteria used 
to define the clusters both in redshift space and in projected separation. 

We have shown (Girardi et al. 1993a, and 1993b) that at least 20 galaxy 
redshifts are generally needed to describe cluster parameters (in particular 
the galaxy velocity dispersion and the virial radius) without heavy problems 
induced by undersampling. So, from our statistical point of view, the choice 



of a minimum of 20 redshift per cluster seems to be an acceptable compromise 
between the need of having at disposal a large number of clusters and the 
possible presence of undersampling problems. 

Girardi et al. (1993a) studied, in a homogeneous way, the velocity field of 
galaxy clusters. They evidenced that the distribution of galaxy velocities in 
clusters is quasi-Gaussian, which encouraged us to apply the Virial Theorem 
to estimate the cluster masses. Moreover, the estimates of cluster masses, 
via the Virial Theorem, have been widely debated and accepted in recent 
literature. As shown by, e.g.. Postman, Geller, & Huchra (1988), Perea, del 
Olmo, & Moles (1990), and Pisani et al. (1992), the Virial Theorem measures 
masses with confidence comparable to that of other methods, in particular 
the Projected Mass estimator (see Heisler, Tremaine, & Bahcall, 1985). 

Our virial masses have been obtained by using robust estimates of the 
velocity dispersions (Beers, Flynn, & Gebhardt 1990, Girardi et al. 1993a). 
The mass errors quoted are obtained via the jackknife method (see, e.g., Efron 
1979, and Geller 1984, for an astrophysical application). In the present paper 
we evaluate the virial masses both using two apertures, 0.75 h^^ Mpc and 1.5 
h"^ Mpc, i.e., half and one Abell radius, respectively. We also computed the 
projected masses of our clusters, within the same aperture. Remarkably, the 



shape of the distributions of virial and projected masses is almost the same. 

The presence of subclustering may bias cluster mass estimates. In par- 
ticular the bias may be severe, and the mass may be overestimated, when 
the subclustering masks the presence of two or more unbound systems. In 
this case one should consider the clumps as isolated clusters or groups. We 
adopted the method of Dressier & Shectman (1988) to investigate the pres- 
ence of subclustering in our sample. Actually this method (hke other methods 
present in the literature) is efficient for a number of members N -^40; for a 
lower number of members the efficiency decreases; therefore we considered 
only well sampled clusters with > 40. The probability of subclustering was 
computed using all the galaxies within 1.5 Mpc from the cluster center; 
a probability > 95 % was taken to be significant. 

The mass distribution of the 18 clusters without evidence for subclus- 
tering did not succeed to be significantly different (according to the Mann- 
Whitney U-test, see, e.g., Siegel 1956) from that of the 12 clusters with evi- 
dence of subclustering. We obtained again the same result when we compared 
the distributions of masses computed inside 0.75 Mpc (27 vs. 16 clusters). 
This result suggests that the effect of subclustering does not strongly affect 
our mass distribution, either because it is small, or because our selection cri- 



teria for cluster membership and our estimate of the velocity dispersions are 
statistically robust (see Girardi et al. 1993a). Therefore, no cluster has been 
rejected. In view of possible surviving cases in which subclustering effects 
produce mass-overestimates, it is conservative to consider mass estimates ob- 
tained within small apertures. As a matter of fact. West & Bothun (1990) 
evidenced that only a few significant subclumps exist within 1 Mpc in 
~ 70 clusters analyzed. 

Of course, other systematic errors in the mass evaluation may be present 
if hght does not trace mass in our clusters (see, e.g., Merritt 1987, Watt et 
al. 1992; and references therein). 

Masses are in solar units. 

2 The Distribution of Cluster Masses 

In the above-mentioned Table 1 we also list the logarithms of the cluster 
virial masses log(/iMo.75), and log(/iMi.5) measured within apertures of 0.75 
h~^Mpc, and 1.5 h~^Mpc. In order to measure mass at a certain radius, 
we must have information on the galaxy velocity-field and the gravitational 
potential up to the same radius. Therefore, we have estimated M0.75 only for 
those 69 clusters which are actually well sampled up to 0.75 Mpc from 



the cluster center; these clusters contain > 20 galaxies all located within 0.75 
h"^ Mpc , and the most external of these galaxies is located, in each cluster, 
at ~ 0.75 Mpcirom the center. Similarly, we have estimated M1.5 only 
for those 47 clusters which are actually well sampled up to 1.5 Mpc from 
the cluster center; these clusters contain > 20 galaxies all located within 1.5 
h"^ Mpc , and the most external of these galaxies is located, in each cluster, 
at ~ 1.5 Mpc from the cluster center. 

The masses M0.75 and Mi. 5 both contribute to describe the cluster mass 
distribution. In fact, fixed apertures may contain different fractions of clus- 
ter masses (depending on the intrinsic cluster sizes), and define two different 
cuts of the density peaks which identify the clusters in the cosmic density 
field. However, in the present paper, we restrict our analysis to the galax- 
ies contained within 0.75 h"^ Mpc , mainly because of the above-mentioned 
problem of subclustering. The distribution of M0.75 is plotted in Fig.l. 

Our sample, although quite extended, is not complete, either in richness 
or in mass. Therefore, we normalized the frequencies of clusters belonging 
to different richness classes to the corresponding frequencies of AGO 's Cat- 
alogue; the normalization procedure is described in Girardi et al. (1993a). 
However, since the completeness of ACO's Catalogue for = clusters is in 



question (see, e.g., Scaramella et al. 1991, and Guzzo et al. 1992), we also 
considered the Edinburgh-Durham Cluster Catalogue (EDCC, hereafter) of 
Lumsden et al. (1992), by re-scaling their richness to ACO's. Moreover, 
using two catalogues also allows us to investigate the sensitivity of the re- 
sults on the choice of the catalogue. In Fig. 2, we show the histograms of the 
ACO-normahzed and EDCC-normahzed mass distributions of our clusters. 

3 The Mass Function 

As one can see in Fig. 2, the mass distributions do not increase monotoni- 
cally from high to low mass values. Two facts may contribute to the de- 
crease at low mass values: (i) the possible incompleteness of ACO's and 
the Edimburgh-Durham Cluster Catalogues for low richness clusters; (ii) the 
lack, in our sample, of rich galaxy groups, the masses of which are expected 
to fall partially in the low- mass region (in fact, completeness in richness does 
not necessarily imply completeness in mass). 

The problem of the non-monotonic behaviour of the mass distribution 
will not be solved unambiguously until a complete sample of galaxy systems, 
spanning the range from small groups to rich clusters, is available. So, even 
if more data on poor cluster and rich groups (which are becoming available) 



were considered, we would still deal with an incomplete sample, and the 
re-normalization problem would not be solved. 

So we deemed it conservative to use only the high-mass (decreasing) part 
of the distribution in order to obtain the mass function of nearby clusters. 
Specifically, we set a lower bound of log(/iMo.75) = 14.6 to the mass distri- 
butions. 

We fitted these two culled distributions with power laws, n{Mi5)dMi5 — 
A(Mi5)"(iMi5, where M15 is the mass expressed in units of lO^^M©. These 
power laws were convolved with Lognormal functions whose dispersions were 
obtained from the median values of the cluster mass errors (quoted in Table 
1). The best-fit values of a (obtained via the Maximum Likelyhood method) 
are — 2.3lo!6 ~'^-^-os for the ACO-normalized and EDCC-normalized 
distributions, respectively. These a values are not different within the errors, 
and the result does not seem to be very sensitive on the choice of the reference 
catalogue. 

The volume density of AGO clusters has been estimated to be 2.5 x 
10~^h^Mpc^'^ for i? > by Scaramella et al. 1991, and Zucca et al. 1993, and 
to be 8.6 X 10-^ h^Mpc-^ and 7.2 x 10'^ h^Mpc'^ for i? = and i? > 1, respec- 
tively, by Peacock & West (1992). Thus the mass function, in the assumed 



range of completeness, is given by n(Mi5) = 5x10 ^/i^-'^(Mi5) "^-^Mpc ^(lO^^M©) ^, 
if one considers the density estimate of Scaramella et al. (1992) and Zucca 
et al. (1993), or ^(Mis) = 3.2 x lQ-^h^-\Mir^)-'^-^Mpc-^{W^MQ)-^ , if one 
considers the density estimate of Peacock & West (1992). 

It is also possible to give an estimate of the exponent n of the power 
law spectrum P{k) oc fc" of cosmic fluctuations. A way to do that is to 
compare the slope of the mass function with that one of the power law repre- 
sentation of the temperature function of clusters given by Henry & Arnaud 
(1991) and Henry (1992). By making use of both of these functions we 
can better constrain the value of n than relying only on the mass function, 
which admittedly span a short range in cluster masses. In the framework of 
self-similar clustering (Kaiser 1990), virial mass M and temperature T are 
linked by T oc M^^"")/^. Comparing the slopes of the temperature function 
(—4.75 ± 0.50) and of the mass functions gives a spectral index n — — l.lioig 
and n — — 1.61J J for AGO and EDCC normalized distributions, respectively. 

In conclusion, we have presented a cluster mass function, obtained directly 
form the dynamics of these galaxy systems, which is described by a power 
law, and which mainly refers to high mass clusters. In the future we plan 
to extend this mass function to a larger mass range, including dynamical 



estimates of masses of poor clusters and groups of galaxies. 
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Captions to Figures 

Fig.l: The cumulative distributions of log(/iMo.75) for our 69 clusters with 
data up to 0.75 h~'^ Mpc . 

Fig.2: The histograms of the mass distributions of our clusters; the upper 
panel shows the ACO-normalized distribution, while the lower panel shows 
the EDCC-normalized distribution. 



